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1. Introduction 

A few years ago, a non-perturbative formulation of Heavy Quark Effective Theory (HQET) 
has been given in [[T]] - see ||^ for a review given at this conference. In particular the problem of 
power divergences is solved through a finite volume matching. Last year, this has been applied to 
the quenched computation of the b-quark mass at the 1/m order In the same spirit, we present 
here a strategy to compute a heavy-light decay constant. We start by writing the Lagrangian at the 
leading order (i.e. in the static approximation) and add a kinetic and a magnetic piece (we follow 
the conventions of [^] and set the counterterm 5„, to zero) 

■i^HQET = Wh ^0 V^h - V^h " (O^^inWh <^ " B V/h • ( 1 ■ 1 ) 

A precise definition of the operators Do, and a-B can be found in Here we just note that 
COkin and Wspin are some bare parameters of the effective theory. 

1.1 Schrodinger functional (SF) correlation functions 

In QCD, we consider the (renormalized and improved) current to boundary correlators /a and 
kv defined - up to improvement factors such as (1 + b/^amf^ \^) - in the SF by 

/a(xo) = -Z^z2^£((Ai)o(x)Cb(y)75Ci(z)), (1-2) 
fcv(xo) = -ZyZ^^ £ ((yi),(x)Cb(y)7^?i(z)), (1.3) 

y,z.k 

where the improved currents A\{x) and V\{x) are defined as in 
We also consider the boundary to boundary correlators 

/i = -Zjf^ £ (C/(u)75Cb(v)?b(y)75?i(z)>, (1.4) 

u,v,y,z 

^i = -2?^ I (C/(u)7.Cb(v)Cb(y)7.?i(z)). (1.5) 

u,v,y,z,<: 

Expanding these correlators at the 1/m order of HQET, and using spin-flavor symmetry, one finds ' 

/a(xo) = zl'^''''z^^Zi^^-'"^'^-^'^[fr{x,) + c''^'^^ , 
^v(xo) = -Z^^''Xz?e"'"{/r(^o)+c^^^V|'in^^ , 

/l = Z|z2e-'"^-^ + O^n/I"" + «spin/r} > 

k, = Z2 z2e-'"-=^ |/f^' + co^jf^ - \co,^^fr^ , 

where mbaie is the (linearly divergent) bare quark mass. 

'The reader can find the definitions of the various correlators /^'f i/ap/aT'-^Ia' ''^ 
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1.2 Basic observable 

We consider a volume Lr' with a time extent T = L, and define (in QCD) 

In the large volume limit, this observable is related to the decay constant, Fb, by 

^■^(D^^^nQFBV^). (1-6) 
In a small volume ^ of space extent L\ ~ 0.4 fm, this observable is matched to its HQET expression 

<I>Q^°(LO=<^^^(LO. (1.7) 

Using the expansions of the correlators /a and f\ given previously, one finds for the rhs at the static 
and at the \/m order 

= InZr + ln +0(l/m), (1.8) 

<j>r(L) wj-{L) = inzr^+in (^^^^) ^'""^""wm ^'-'^ 

( f^L/i) ff^\ ( fr{L/2) fr\ ,^.. ,2. 

1.3 Evolution to larger volumes, in the static approximation 

In order to clarify the discussion, we first explain the strategy in the static approximation, the 
generalization to the \/m order will be done in the next section. We start by the matching of HQET 
to QCD in the volume Li, at the static order : <I>p™(Li) = <I>p^'(Li). The evolution to a volume 
Loc = L3 = 2L2 = 4Li is then done, within the effective theory, in the following way: 

<I>F(Lec) = [^f\L^)-^f\L2)] + [^f\L2)-^f\U)]+^f''{U) . (1.10) 

We note that Z^^^ cancels in the differences [<I>p^'(2L) — <I>p^'(L)]. Using the renormahzed SF 
coupling g^{L) [Q], we define the static step scaling function (ssf) 

ar{u) = [<I>r(2L)-<I>r(L)]_,(^^^^ = Jim^ir(",a/^) (1-11) 



Er(^,a/L)= [r'(2L)-r'(L)]_,(,^^„ where r'(L)= In ^^11 ^ . (1.12) 



-friL/2)\ 

We can now rewrite the rhs of eq ( |1.10| ) as the sum of three continuum terms 

^piL^) = ar(M2) + af '("i) , where = fiU) • (1-13) 

Before discussing the 1/m corrections we close this section by a few remarks : 

^The matching is done in a small volume, in order to be able to simulate a b-quark with the discretization errors 
under control. 
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• The different terms in the eq ( 1. 12) have to be computed at the same value of the lattice spac- 
ing (in order to insure that has a well defined continuum limit). This is due to divergences 
proportional to the logarithm of the lattice spacing that one has to cancel. 

• In eq. ( 1.13 ), the entire quark mass dependence comes from <I>p^^(Li). 

• At this order, since there is only one matching constant to eliminate (Z^^'), it is sufficient to 
match one observable (<I>f)- 

1.4 Including 1 /m corrections 

At this order, there are three more matching parameters in Op compared to the static case ^ . There- 
fore, to determine them, we introduce three other observables defined in a volume of space extent 
L and we give their expressions at the \/m order 

(L) ^\{R\ + 3/?7) - = OhinRf"" , (1.14) 

^2{L) = -ln ij^J = cospin^ with T=L/2, (1.15) 

<^,{L) ^ RA{L/2)-Rf\L/2) = c«Q^Xa(L/2) + (Oki„/?^"(L/2) + «,pi„/?f "(L/2) ,(1.16) 

where the definitions of the ratios R can be found in ^ ^. Together with Op, given at this order 
by (L9), we then have a set of four observables. In these observables, we have chosen to subtract 
the static part (when existing) from the QCD one, as we did in the mentioned reference. This is 
perfectly legitimate because they both have a continuum limit, and this simplifies the equations. 
Like in the static case, the matching is imposed in the volume L\. This allows us to replace in (1.9) 
the parameters (Okin, (Ospin and c^^^^ by a combination of QCD and HQET quantities. 
The evolution to the volume L2 is given by 

^AL2) = <I>[?Q^'^(L2)-<I>|?Q^'^(L0 + <I>r(Li) (1.17) 

= [<I>r(L2) -OfHZ-i)] + [^'J'"{L2)-^'/"{L,)] +<I>f °(Li) . (1.18) 

The ssf for the static term has already been given in the previous part, and for the l/m part we write 

^'J''\2L)-^'/\L)=l^a,{f{LmiL) . (1.19) 

i 

The expressions for the ssf can be found from the last equation by using ( |1.9| ) together with (1.14), 

( 1.15 ), ( 1.16 ) in the volume Li. The explicit definitions are given in the appendix. 

In the step L2 — > Loo, we need <I>i(L2,M), and we are then lead to define the ssf for the <I>i : 

«I>i(L2)=£aij(f2(L))<I>j(Li). (1.20) 

We can write down the final equation for Op 
<Dp(L„) = ar("2) + af'(Mi) + XcTi(«2)aij(Mi)<I>j(Li) + £cTi(Mi)<I>i(Li)+<I>Q^°(Li) (1.21) 



3 Also Z^^^'^ is different thaii Zf\ but as in the static case, it drops out in the differences 

^We remind the reader that the quantities defined with an subscript 1 are "boundary to boundary" observables. 
Because the noise over signal ratio grows exponentially with the time, we impose T = L/2 for all these observables. 
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2. Results in the quenched approximation 

We used basically the same data as in in which the reader can find the details of the 
simulation. The simulations are done with non-perturbatively 0(a) -improved Wilson fermions. 
The light quark mass is fixed to the strange one. Concerning HQET, we used the HYP actions [^], 
which help to have a reasonable statistical precision. We show the continuum extrapolations in 

fig-in 

The extrapolations are done linearly in (a/L)^ for QCD as well as for the static part, but ina/L 



o 0.6 
0.4 
0.2 

0- 



15 

X 10"' 




Figure 1: Continuum extrapolations of the various terms appearing in eq. ( 1.21 ). The QCD contribution 
Op*^'^ (diamonds) and the static part are shown on the left (the circles represent the large volume contribution 
Of'^'(m2), and the squares the small volume part o^''{u\) ). The 1 /m correction is plotted on the right, the 
circles represent the large volume terms Lij 0'i(M2)0ij(Mi)^'j(ii), and the squai^es the small volume part 
Ii(7i(Mi)<I>i(Li). The color blue stands forHYP2, and red for HYP 1 [||. 

for the \/m term, because of the absence of 0(a) improvement. The ordinate scale is the same in 
order to compare the relative size of the different contributions. Concerning the precision, one can 
see that the total error is largely dominated by the one of the \/m part in the large volume. Since 
for this part the results are not yet completely satisfactory, we refrain from performing a continuum 
extrapolation. We will use the result at the finest lattice spacing only (j3 ~ 6.45, L/a = 32). 
Our preliminary results for Fb^ are shown in table [T[ In the first column we give the results in 
the static approximation, while in the other columns, we have included the \/m corrections. We 
observe that in the static approximation, depending on the matching condition represented here by 
6 ^ , the result can change by 7%. This variation disappears when the l/m terms are included. 
Note that differences of f_^^^+^/™^ table |l|, have much smaller error than their individual values, for 
example 

^^tat+l/m^g^ = 0, 01 = 1, 02 = 0) -Fg "'+'''"'(00 = 1, 01 = 0, 02 = 0.5) = 4±2 McV. (2.1) 

The other information is that the l/m term contributes (with a minus sign) up to ~ 15% to the 
final result. One can see that adding the l/m terms increases the size of the statistical errors, as 
expected from the previous plots. This is due to the fact that the signal for the l/m part in large 



The quark fields are periodic in space up to a phase 9. 
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volume is more difficult to extract than in the static case, and also because of the absence of 0(a)- 
improvement at this order. We also note that our result is compatible with a recent computation 
done with a different method, but which also goes beyond the leading order of HQET [^. 





rstat 


j-,stat+l/m 


00 




01 =0 

02 = 0.5 


01 =0.5 

02 = 1 


01 = 1 

02=0 





224 ±3 


185±21 


186 ±22 


189 ±22 


0.5 


220 ±3 


185±21 


187 ±22 


189 ±22 


1 


209 ±3 


184±21 


185±21 


188±22 



Table 1: Results for Fb^ in MeV with and without the 1/m corrections, for different values of the 9 angles. 



3. Conclusion 



We have shown how to perform a non-perturbative computation of a heavy-light decay con- 
stant at 1/m order of HQET, and we have given preliminary numerical results in the quenched 
approximation. The inclusion of the dynamical quarks is on the way [^, ||]. Applying this method 
for Nf > should allow for precise computations of the heavy-light decay constant, with a good 
control on the systematic errors. Note in particular that eq. is a good sign of the absence of 
significant 1 /m^ corrections. On the numerical side, the cancellations of the divergences require 
sufficient statistical precision, and we hope that the all-to-all propagator, like proposed in [^, [l^ 
will be of great help there. 
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4. Appendix: The step scaling functions 

In order to have more compact notations, we replace Op by O4, such that the ssf Oi introduced 
in eq. ( |1.19| ) are now represented by 041 . We can rewrite eq. ( |1.19D together with eq.( |1.20D as 



<I>i(Lk+i,M) = ^ aij(Mk)<I>j(Lk,M)±54ia^(Mk). 



(4.1) 



The ssf are then given by a four by four matrix 
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022 










032 







041 


042 


043 


1 
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To give their explicit expressions, we define 



I ) V/|tat(r/2) fl 

Then, one finds: 



) = \pr{2L)/pr\L: 



and ps, 



I3i(m) 



S32(«) 



£41 (m) 



£42 (m) 



£43 (m) 



1 fj,^n.,r ^ _ /?^"(^)^5a(2Z.) 

1 



priL) 
1 



pr(L) 



f(L)=u 



RsAL) 



f{L)=u 



»I/kin(2L)-»F'^"(L) 



»I/spin^2L)-»F'P'"(L) 



PS^{2L) - ps^{L 



Rf{L) 



(kin , X iT.kin \ PSa (2^) ~ PSa i^J nkin 



^5aW 



f(L)=u 



PS^{2L) - ps^{L) 



RsAL) 
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